Probabilities in the landscape: The decay of nearly flat space 
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We discuss aspects of the problem of assigning probabilities in eternal inflation. In particular, 
we investigate a recent suggestion that the lowest energy de Sitter vacuum in the landscape is 
effectively stable. The associated proposal for probabilities would relegate lower energy vacua to 
unlikely excursions of a high entropy system. We note that it would also imply that the string theory 
landscape is experimentally ruled out. However, we extensively analyze the structure of the space 
of Coleman-De Luccia solutions, and we present analytic arguments, as well as numerical evidence, 
that the decay rate varies continuously as the false vacuum energy goes through zero. Hence, low- 
energy de Sitter vacua do not become anomalously stable; negative and zero cosmological constant 
regions cannot be neglected. 



I. INTRODUCTION 

In any theory with more than one metastable or stable 
vacuum, the question arises with what probability cer- 
tain low energy physics phenomena will be observed. It 
is not enough to count vacua Q, |^ |1| with the speci- 
fied property and fold in anthropic constraints, because 
the cosmological evolution may favor some vacua over 
others. Whether and how inflationary expansion factors, 
decay rates, and initial conditions affect the probability 
to be in a given vacuum is a major challenge in theoret- 
ical cosmology. Even once st ring theory and its vacuum 
strurturejO 11 11^^ ip, [H El 111 Eljll El 
ElElElilillllilil, or "landscape" Hg, are 
fully understood, the theory will become predictive only 
if this problem is solved. 

At the core of the problem is the phenomenon of eter- 
nal inflation. In a global description, an infinite num- 
ber of infinitely large regions containing each vacuum 
are produced, after unlimited volume expansion of the 
preceding vacua. Elegant proposals to regulate these in- 
finities have been advanced (see Ref. |2g for a recent 
review, and Ref. |27j | for an early discussion), but so far 
no prescription stands out uniquely, and each has some 
counterintuitive properties^. 
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^ For example, recent proposals |2SL l2fll . when applied to a land- 
scape consisting of two de Sitter vacua, would rate them equally 
likely even if their cosmological constants (and thus their en- 



Proposals for probabilities, and indeed the phe- 
nomenon of eternal inflation itself, are usually formulated 
in terms of a global, scmiclassical geometry. In a separate 
publication |3l| , it will be argued that such a description 
does not exist. Independently of this point, however, one 
has reason to be sceptical of any simultaneous descrip- 
tion of regions that are forever causally separated. In the 
context of black hole formation and evaporation, it is im- 
possible to reconcile a global description with quantum 
mechanical unitarity and linearity |82l |. A good theory 
should be capable of describing any one observer's ex- 
perience, for example an observer inside the black hole, 
or an observer remaining forever outside. A description 
of two causally disconnected observers at once, however, 
does not correspond to any feasible experiment, and if 
attempted, leads to contradictions. Adaptations of this 
lesson to the cosm olog ical horizon of de Sitter space in- 
clude Refs. Hi, mill El El El SHI. 

This suggests that eternal inflation should be reformu- 
lated in terms of a single causally connected region, or 
causal diamond [3^ . (A causal diamond is the overlap of 
the causal past of q with the causal future of p, where p 
and q are any two points on a worldline. The maximal 
area on the boundary of a causal diamond is an upper 
bound on the entropy in the region it encloses |3H40l|.') 
Perhaps a compelling deflnition of probabilities can be 
given in this observer-centered language. 



tropy) differ by an enormous factor. In an older proposal l3Cl , all 
metastable vacua (including, presumably, our ow n) w ould have 
vanishing probability if any stable vacuum exists \2f\ . 
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This program is straightforward to implement for a po- 
tential landscape whose lowest point has positive vacuum 
energy. Because the entire system has finite entropy, it 
will probe its whole phase space over and over, and the 
relative probabilities are easy to obtain. Under unitary 
ergodic evolution, the relative amount of time spent in 
a given vacuum, and thus its relative probability, is pro- 
portional to the number of states it contains, Ni. (We 
neglect anthropic factors in this discussion.) 

The entropy a de Sitter vacuum is given by the area of 
its cosmological horizon: iS = Ai/A, where Ai = 12n / h.i. 
Hence, 

M = exp(37r/Aj) (1) 

for a vacuum with cosmological constant A^. The lowest 
vacuum has the largest entropy and can be thought of 
as a highly degenerate ground state for the landscape. 
Its exitations correspond to less entropic configurations, 
such as a nonempty de Sitter space. The other vacua are 
also merely particular excitations of the lowest vacuum. 

While a landscape with only positive energy vacua is 
pleasantly simple, it is also experimentally ruled out. 
Dyson, Kleban, and Susskind |32j showed on statistical 
grounds that our universe cannot have arisen from the 
unitary, ergodic evolution of a completely stable de Sitter 
spacetime. Their analysis can be adapted to any land- 
scape whose lowest point has positive vacuum energy, 
with the same conclusion. Unless one of the assumptions 
(unitarity, ergodicity) breaks down, this implies that the 
true landscape contains vacua with nonpositive cosmo- 
logical constant. 

Indeed, the string theory landscape is expected to con- 
tain valleys with negative cosmological constant as well 
as supersymmetric regions with vanishing vacuum en- 
ergy jj]. The latter pose a significant challenge for a 
causal definition of probabilites, in that they tend to at- 
tract all of the probability, at least according to very 
simple probability measures one might propose.^ (This 
is reminiscent of the difficulties with the global prescrip- 
tion of Ref. ^3 , where the probability has support only 
on nonpositive vacua.) 

To see this, consider two very simple proposals. First, 
let us suppose that the probability for a given vacuum is 
still proportional to the amount of time one can spend in 
it. This probability measure has support only for vanish- 
ing cosmological constant. This is because all the de Sit- 
ter vacua decay after a finite time, and all the vacua with 
negative cosmological constant suffer collapse into a big 
crunch after a finite time of order |Ai|~^/^. The A = re- 
gions, however, are open FRW universes that live forever, 
so they dominate. 

Second, suppose that the probability of each vacuum is 
set by the maximum entropy. More precisely, suppose it 



is proportional to the number of quantum states that can 
exist within a causal diamond contained in a correspond- 
ing region. The maximum entropy is set by the largest 
area on the boundary of the causal diamond , which 
goes like an inverse power of the cosmological constant. 
In particular, one finds that this area is finite for de Sit- 
ter vacua (because of the cosmological horizon), as well 
as for negative energy vacua (because of their finite life- 
time). On the other hand, the causal diamonds in open 
FRW regions with vanishing cosmological constant can 
become arbitrarily large, so they dominate. 

This would seem to suggest that a more refined proba- 
bility measure is required. Banks and Johnson |46l | have 
recently proposed a different approach, in which the de- 
scription of the full string theory landscape would be 
quite similar to that of the above toy landscape with 
only positive energy vacua. The lowest positive energy 
vacuum would play the role of a degenerate ground state 
in which the system spends most of its time. ^ 

In support of this proposal. Banks and Johnson pre- 
sented numerical and analytic arguments suggesting that 
very low energy de Sitter vacua are anomalously long- 
lived. That is, a vacuum with Agmaii much smaller than 
the barrier height would be significantly more stable than 
a zero-cosmological constant vacuum obtained by shift- 
ing the entire potential down by the tiny amount Agmaii- 
More precisely, the lifetime of the Agmaii vacuum would 
diverge as Agmaii approaches zero from above, but would 
become finite again when Agmaii = 0. 

We feel that the Banks- Johnson proposal does not re- 
move the difficulties with A = vacua, whose entropy 
and lifetime can be arbitrarily large. Moreover, the pro- 
posal involves ad hoc assumptions that contradict semi- 
classical results in a regime where the latter are expected 
to be valid. This includes the unavailability of vacua with 
small negative cosmological constant, and the ability of 
negative energy vacua to "decay back up" to the putative 
de Sitter ground state so as to establish detailed balance. 

But if these obstacles could be overcome, the Banks- 
Johnson proposal would in fact rule out the string theory 
landscape. (This was not noted in Ref. '4^.) This is be- 
cause the arguments of Dyson, Kleban, and Susskind [39| 
would then apply, and the universe we observe would be 
overwhelmingly unlikely to arise dynamically. 

In any case, the suggestion that low energy de Sitter 
vacua are anomalously stable is of independent interest 
and warrants careful investigation. It is a radical claim 
in that it conflicts with locality. When a bubble of true 
vacuum appears it is typically of microphysical size, and 
one would not expect it to matter whether the universe 
has a cosmological horizon out at, say, 10'' Mpc. Yet, 
according to Ref. 46], it matters a lot: the addition of a 
small amount of vacuum energy (the smaller, the better!) 



^ For the sam e reason, they also promise to admit precise physical 
obseryables I4ll4ll4ll43 . 



^ Another discussion of the Banks-Johnson proposal can be found 
in Ref. 
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would drastically suppress the decay of flat space. ^ 

In this paper we find that the decay rate is contin- 
uous as the false vacuum energy passes through zero.^ 
Hence the low energy vacua are not anomalously long 
lived, i.e., not much longer lived than the flat space vacua 
they would be shifted into. 

Depending on the detailed shape of the potential, a 
A = vacuum can be completely stable. Then conti- 
nuity demands that the decay of a false de Sitter vac- 
uum related by a shift of the potential become arbitrarily 
suppressed as the false vacuum energy app roaches zero. 
(The numerical data presented in Ref. |4a| actually per- 
tain to this case and thus did not support the conclu- 
sions drawn there.) Thus, our results do not exclude the 
possibility of extremely long-lived de Sitter vacua (up to 
log^docay ^ What we rule out is that they will 

generically be far more stable than the flat space vacua 
obtained by shifting the potential down. 

It may be that the A = vacua in the landscape are 
necessarily supersymmetric. (Why else would the contri- 
butions to the cosmological constant cancel precisely?) 
Then they would all be stable. However, this does not 
automatically imply that the low energy de Sitter vacua 
have lifetimes comparable to the recurrence time, since 
it is not clear in which sense those vacua are close to su- 
persymmetric. In the real landscape, we do not get to 
shift the potential continuously. A low energy de Sitter 
vacuum may be very far from a supersymmetric region 
in fleld space. 

The observation of Banks and Johnson stands that the 
decay of de Sitter space is exponentially suppressed while 
negative energy regions meet their demise in polynomial 
time. The potential significance of this asymmetry re- 
mains to be explored. It would be premature to con- 
clude, however, that negative energy vacua can simply 
be neglected. In any case, a more immediate challenge 
is that the A = vacua appear to dominate, both in 
terms of entropy and in terms of their lifetime. This is 
problematic because we do not appear to live in one. 

This paper is structured as follows. In Sec. we re- 



Banks and Johnson suggested a holographic explanation of the 
alleged discontinuity. For small positive A, the false vacuum has 
enormous entropy due to the cosmological horizon. The argu- 
ment is that the decay is entropically suppressed, because all 
the horizon entropy would be destroyed by the decay. No such 
contribution to the entropy is present at A = 0, so there should 
be no entropic suppression of the transition to a lower vacuum. 
However, it is not clear to us how adding entropy should stabi- 
lize a system. The cosmological horizon is far away when the 
bubble first forms, and will only be destroyed once the bubble 
has expanded far enough. Indeed, the same logic would imply 
that the decay of fiat space could be suppressed by adding enor- 
mous entropy to faraway regions (which can be done by adding 
a huge black hole, or with arbitrarily small backreaction using 
radiation). 

^ This was discovered independently by A. Aguirre, T. Banks, and 
M. Johnson, whose results appear simultaneously with this pa- 
per. 



view the instanton calculus describing the decay of a false 
vacuum in the presence of gravity. In Sec. lIIII we demon- 
strate the continuity of the decay rate. The key ingredi- 
ent in our argument is to show and exploit that singular 
solutions nearby the regular instantons behave continu- 
ously as the false vacuum energy is taken to zero, unlike 
the regular instanton geometry itself. 

In Sec. II VI we give a thorough analysis of the space of 
solutions to the Coleman-De Luccia equations, depending 
on the initial value of the field, and on the energy of the 
false vacuum. We now also explore regions where the 
latter is far from zero. In Sec. we present numerical 
evidence supporting our analytical results. 



II. COLEMAN-DE LUCCIA TUNNELING 

In this section, we briefly review the decay of a false 
vacuum, following the analysis of Coleman and De Luc- 
cia 01 . (For earlier work that neglects effects of gravity, 
see Refs. ^3^|E9-) We set up notation, mostly following 
Banks and Johnson 46]. 

Consider a scalar field (j> with potential V{4>). We as- 
sume that V has two local minima, at (jjp and with 
V{(I)f) > V{(j)T). Examples are shown in Fig. 01and Fig. El 
below. Without loss of generality, we take the top of the 
barrier to be at </) = and the true vacuum at 0t > 0. 

Let us assume that initially the field is in the false 
vacuum thoughout space: (f) = (j)-p. Classically, it would 
remain there forever. Quantum mechanically, it may be 
possible for it to lower its energy by tunneling through 
the barrier towards the true vacuum. The path integral 
for this process can be approximated by a regular Eu- 
clidean solution, or instanton. 

The most symmetric (and, presumably, dominant) 
nontrivial Euclidean solution is S'0(4) invariant, with 
metric 



ds^ = dt^ + p{tfdn 



3 ' 



(2) 



Here dJlg is the metric on the unit three-sphere; we will 
think of the parameter t as Euclidean time. Thus the 
instanton is described by two functions 0(i) and p{t). 

The equations of motion are Euclidean versions of the 
FRW equations for a closed universe with a scalar field: 



P - 



An 



0' + vi4>) 



p 



(3) 
(4) 



where an overdot (prime) denotes differentiation with re- 
spect to t {(/)). They obey the constraint 



Stt 



(5) 



Thus, the scalar field behaves like a particle moving in the 
Euclidean potential U = —V (see Fig. QJ, with friction 
proportional to p/ p. 
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There are compact and noncompact solutions. Non- 
compact solutions have one pole (a value of t, conven- 
tionally taken to be zero, at which p vanishes), and have 
the topology of R^. Compact instantons have two poles 
and have the topology of a four-sphere. For compact 
solutions, p, and thus the friction, eventually becomes 
negative. 

At any pole, continuity of the field gradient requires 



{at p = 0) 



(6) 



Moreover, the absence of a conical singularity requires 
the boundary condition 



\p\ = 1 (at p = 0) 



(7) 



At the first pole, from which the equations are integrated, 
this is imposed as a boundary condition. For compact 
solutions, the constraint equation ((SJ guarantees that 
Eq. will automatically be satisfied at the second pole, 
as long as it has been ensured that the scalar field energy 



(8) 



remains bounded. For compact solutions, Eq. |(7J) implies 
divergent antifriction at the far pole. 

Let us first consider the case where the false vacuum 
has zero or negative energy: V{(j)F) < 0. Decay is me- 
diated by an instanton that is asymptotic to the back- 
ground, i.e., to Euclidean flat or AdS space. Hence, the 
instanton must be noncompact, with — > (/)f for p — > oo. 
It must contain a bubble of true vacuum: > at p = 0. 
As we shall discuss in more detail below, it is by no means 
automatic that such an instanton exists. If it does not, 
the false vacuum is stable. If it does, then the rate of 
decay per unit volume is given by 



r exp[-(5inst - Sh 



(9) 



Here, S'inst is the Euclidean action of the instanton, and 
^bg is the action of the background solution. (This is 
the solution corresponding to 0(t) — 0f for all t, i.e., 
Euclidean flat or AdS space.) 
The action is given by 



S = d Xy/g 



PI 



IBtt 



R 



1 



(V0)2 + y(0) 



On shell, this becomes simply 



The energy obeys 



E 



(10) 



(11) 



(12) 



Since p > for noncompact solutions, it is clear that 
energy will be lost to friction. Hence, the scalar must 



start out, at p = 0, at some value 0o where the Euclidean 
potential is higher than in the true vacuum: U{(j)Q) > 
U{<j)T). A suitable point may or may not exist. By the 
same token, it is clear that there is never an instanton 
describing the formation of bubbles of false vacuum in a 
true vacuum with nonpositive energy. 

Now let us turn to the case where the false vacuum is 
de Sitter space: T^(0f) > 0. (The above equations for the 
action and the energy still stand.) In this case the back- 
ground action in Eq. is finite, and tunneling is always 
allowed. (If all else fails, the Hawking-Moss instanton, 
0=0, describes tunnehng to the top of the barrier [5H| : 
see also Ref. |E2-) However, if the background action 
dominates over the instanton action, tunneling will be 
extremely suppressed: logF ~ exp5bg, corresponding to 
a lifetime of order the Poincare recurrence time in the 
background de Sitter space. 

The instanton, in this case, will be compact and will 
not reach either 0t or 0f exactly. As we will discuss in 
more detail in Sec. IIVI there can be several such instan- 
tons. We will be most interested in one-pass instantons, 
which cross the barrier precisely once. An example (type 
I) is an instanton that resembles the de Sitter four-sphere 
background solution except in a small region containing a 
bubble of true vacuum. Then the action difference results 
from the bubble region alone, and the rate will not de- 
pend strongly on the background cosmological constant. 
Another example (type II) is an instanton that does not 
spend much time near (p-p. Then the background action 
dominates, and tunneling is suppressed by the inverse 
background cosmological constant in the exponent.^ 

In the following we will investigate how the tunneling 
rate depends on overall shifts of the potential. We write 



v{4>) = v^{<t>) + Vp , 



(13) 



where Vo(0f) = 0. The additive constant V-p shifts the 
entire potential up or down, with Vp being the false vac- 
uum energy. 

Although we will explore the parameter space quite 
broadly, we are interested mainly in confirming the con- 
tinuity of the rate as the false vacuum energy passes 
through zero. For potentials V that permit the decay 
of flat space, we will find that the addition of a small 
cosmological constant V-p does not change the rate much, 
leading to a type I decay. If flat space is stable, then 
the addition of a small cosmological constant leads to a 



® If the both vacua have positive energy, then tunnehng is possible 
in both directions, mediated by the same instanton. The differ- 
ence in rates comes entirely from the different background actions 
that must be subtracted depending on the direction of the pro- 
cess. From Eq. ^ one finds that Ti^/Ty ~ exp(5'T — 5'f). The 
action of Euclidan de Sitter space is minus the entropy {S = —S). 
Thus, this result agrees with statistical expectations, as discussed 
in the introduction. Unless the two vacuum energies are very 
similar, the upward decay is an example of a highly suppressed 
process where the background action dominates. 
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highly suppressed, type II decay of the resulting de Sitter 
space. These results express the continuity of the decay 
rate. 

In the next section we develop analytical arguments 
in support of these assertions. Global properties of the 
solution space will be discussed in Sec. IIVI The final 
section contains numerical evidence. 



III. CONTINUITY OF THE DECAY RATE OF 
NEARLY FLAT SPACE 

In this section, we present our argument that the de- 
cay rate is continuous as the potential is shifted by a 
constant and the false vacuum energy approaches zero 
from above {Vp = 0). Our argument combines a num- 
ber of intermediate results, which will also be useful in 
later sections. Hence we will structure this section by 
developing individual results in separate subsections. 

We will rely heavily on the properties of singular solu- 
tions. Hence, from here on, when we say "solutions", this 
includes singular solutions, i.e., integrals of the equations 
of motion that run into a singularity. We will write ex- 
plicitly "singular solutions" or "regular solutions" to re- 
fer specifically to one of these classes. "Instanton" means 
"regular solution". When we say "compact" , we mean 
that the radius approaches a second zero, independently of 
whether it does so in a singular or regular way. Through- 
out the paper, we use the term "generic" in the technical 
sense: namely, a feature is generic if it occurs in an open 
set of moduli space. In other words, generic means "un- 
changed by infinitesimal perturbations. " 



A. Solutions are generically compact, and all 
singular solutions are compact 

It is important to our method that solutions are gener- 
ically compact, so we will give a careful argument. By 
construction, all of our solutions start with regular initial 
conditions at p = 0; then we evolve in Euclidean time to 
find the full solution. We claim that generically p returns 
to 0; in particular, it will do so for all singular solutions. 

For the sake of this argument, we assume that outside 
the region of interest the Euclidean potential is negative 
everywhere and has no extrema, and also that it stays 
finite everywhere. Recall that we are assuming that the 
local minimum in the Euclidean potential is negative, and 
that the potential has no other extrema. 

There are three ways one might imagine that the radius 
p can fail to return to zero: 

(1) p could remain finite and the solution could remain 
smooth as t — > oo. 

(2) Evolution could stop at a singularity at finite p. 

(3) p could asymptotically approach infinity as i — > cx). 

(4) p could diverge at finite t. 

We will see that possibilities (1), (2), and (4) are for- 
bidden by the equations of motion, while possibility (3) 



occurs nongenerically. An example of option (3) is the 
usual instanton which mediates the decay of a false non- 
de Sitter vacuum. 

To see that possibility (1) is forbidden, we focus on the 
evolution of the Hubble parameter, H — pj p. Combining 
the FRW equations © gives 



H = -47rM, 



-2 12 



PI 



(14) 



The right side is negative semidefinite, so H can only 
decrease. 

We argue by contradiction. Suppose there exists a 
maximum radius Pmax- Then H is bounded above by 
^1/Pinax- This implies that H ^ — oo ast— s-cx). Now 



dt'p{t') 



dt'p{t')H{t') 



(15) 



By assumption, p{t) stays finite, and it must be positive 
because it is the radius of the 3-spheres. Since we have 
just shown that H diverges to negative infinity and by as- 
sumption p{t) is finite and positive, the integral diverges. 
This contradicts our assumption that p remains finite. 

Possibility (2) is a singularity at finite p. To achieve 
a singularity at finite p, we would need p and/or (p to 
diverge. We choose the singularity to occur at i = and 
characterize the divergent behavior of p and (f) there. In 
order for p to diverge without p diverging, we need p ^ t^ 
for > p > -1. The FRW equation is 



1 = 



Stt 



(16) 



By assumption, V is finite everywhere and the singularity 
forms at finite p. Keeping just the singular terms the 
equation becomes 



Stt 
6M2j 



2 12 



(17) 



implying that diverges in the same way as p, ~ t^. 
Now the equation of motion for cj), keeping divergent 
terms, is 



3^, 
P 







(18) 



With our assumptions so far, the first term diverges as 
i^"^, while the second term diverges as t^'^ . Since we 
require > p > — 1, the two terms diverge as different 
powers of t and the equation cannot be solved. 
Possibility (4) is forbidden because 



p(t) = p(to) exp ( / H{t')dt' 

'to 



(19) 



By Eq. (|14|) . H is monotonically decreasing. Thus, if p 
and H are finite at some time t^, then p cannot diverge 
in finite time. 
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This leaves only possibility (3), a regular solution in 
which p —f oo. This is not forbidden, but as we will now 
show, it is nongeneric, as it requires infinite fine-tuning 
of the initial conditions. The basic idea is that if the ra- 
dius p is to approach infinity, then the field will have to 
asymptotically approach one of the vacua so as to pro- 
vide a resting place with negative cosniological constant. 
Since the field must come to rest at a local maximum of 
the Euclidean potential, this requires infinite tuning. 

To show that the field must approach a vacuum, let us 
recall that H is negative semidefinite and H must remain 
nonnegative for a noncompact solution. Hence, H must 
approach zero as p —> oo. Further, Eq. 114() indicates 
that H ^ requires — > 0. Using this information, the 
Friedmann equation Eq. ISJ becomes, as p —> oo, 



Svr 
3M2j 



1^(0), 



(20) 



which requires the potential V to be negative. 

Finally, recall that the equation of motion for (f> is 



(21) 



We have already established that 0, so we can ignore 
the second term. In order for cj) to remain zero, we need 
(j) 0, so V'{(j)) 0. Hence the field approaches a 
critical point of the potential. Since we also showed that 
F < as p — > cx), this means that asymptotes to one 
of the vacua. (As discussed in the next section, this is 
necessarily the false vacuum, though this is not crucial 
here.) 

Now it is easy to argue that this solution is nongeneric. 
Approaching the top of a hill as p ^ cxd is highly unstable. 
An infinitesimal perturbation of the starting point will 
cause the field to overshoot or undershoot the top of the 
hill, leading to a compact singular solution. 

This completes our argument that compact solutions 
are generic. We have also seen that noncompact solutions 
are necessarily regular. Hence, all singular solutions are 
compact. 



B. Noncompact regular solutions are one-pass 

Noncompact solutions have no anti-friction; as shown 
above, H > throughout the solution. Since there is no 
anti-friction, the energy can only decrease. If the field 
starts near the false vacuum, it has no hope of achieving 
the true vacuum, which is at a higher Euclidean energy. 
If it starts near the true vacuum, it cannot have a turning 
point because a turning point requires E < Up < Ut, so 
the field will never have enough energy to ascend to the 
top of either vacuum. 

To summarize, the only possible noncompact solution 
is a one-pass, regular solution with the field near the 
true vacuum at the origin and asymptoting to the true 
vacuum at infinity. 



For oddly shaped potentials, there may even be more 
than one such solution. Such potentials can retain their 
odd properties under small deformations, so using our 
technical definition we cannot call them nongeneric. 



C. The field escapes to ±oo at the singularity 

We showed above that all singular solutions are com- 
pact, with the singularity developing as the second zero 
of p is approached. It follows that the behavior of the 
fields near the singularity is universal. In particular, the 
field (j) escapes in the direction indicated by the values of 
(j) sufficiently close to the second pole. 

To show this, recall that the potential is finite every- 
where. We write the equations of motion, keeping only 
terms which will diverge near the pole. 



3H(P 



PI 







(22) 
(23) 



To find the divergent behavior, we assume that (j) and p 
are both power-law in t. We find 



i{t) = At^/^ 4,{t) = B/t 



(24) 



where we have taken the pole to be at t = 0. Note that (j) 
remains of constant sign, as advertised. The field value 
(j) diverges logarithmically. 



D. Across a regular compact solution, the number 
of passes generically changes by one. 

Consider perturbing the starting point of a compact 
instanton. If the perturbation is very small, then the 
new solution is virtually unchanged for a long time. But 
we know that near the opposite pole the field must blow 
up. This is because the new solution will fail to approach 
the pole with just the right conditions to achieve (j) — 
at the pole, and diverging anti-friction will magnify the 
mistake. 

Generically, perturbing the starting point in opposite 
directions will result in opposite runaway of the field at 
the far pole. A perturbation in one direction will re- 
sult in 4> reaching zero already for finite (small) p. Then 
the field changes direction, and the singular behavior of 
Eq. (|24|) will push it off to infinity. Hence, it will have 
one additional pass compared to the regular solution. 

Under a perturbation of the starting point in the other 
direction, cj) will not reach the zero that it reached for 
the unperturbed instanton at the far pole. The diver- 
gence will push it off to infinity with no additional pass 
compared to the regular solution. Hence, the number of 
passes changes by one across a regular compact solutions. 

Things are more complicated for noncompact instan- 
tons. Because there is no anti-friction, we cannot appeal 
to the universal divergent behavior of Eq. H24|l . Hence, 
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the details of the potential can be important. There is no 
general rule about how the number of oscillations changes 
across a noncompact instanton. Indeed, we will provide 
analytic and numerical evidence that jumps can be by 
more than one oscillation. 

Going in the other direction, we expect that if two 
nearby starting points 0i and 02 result in singular so- 
lutions with different number of passes pi and p2, there 
will be at least one regular solution with starting point 
between (/>i and (fi2- However, we do not prove this rig- 
orously. 



E. The decay rate is continuous near Vp = 

What makes it hard to rule out a discontinuity in the 
tunneling rate at = is the fact that the regular 
instanton does change form: for Vp > it is compact, 
while for = it is noncompact. As a result, doing 
perturbation theory in Vf around the point Vf = is 
confusing. We circumvent this problem by perturbing 
the singular solutions near the instanton. 

There are two cases. We defer the case where flat space 
is stable, and begin with the more interesting case where 
it can decay. Thus, we assume that there exists a non- 
compact instanton mediating the decay at Vp = 0. We 
want to prove the existence of an instanton for infinites- 
imal positive whose action scales with the de Sitter 
radius in the false vacuum. 

Consider perturbing the starting point of the Vp = 
instanton by a small amount S(f>. The instanton itself is 
a solution in which the field starts near the true vacuum 
at yO = and approaches the false vacuum as p oo. 
Generically, perturbing the starting point in one direc- 
tion will cause the field to overshoot the false vacuum, 
leading to a singular one-pass solution, while perturbing 
the starting point in the other direction will cause the 
field to undershoot the false vacuum, resulting in a singu- 
lar solution with at least two passes.^ The instanton is at 
the boundary between the singular one-pass overshoot- 
ing solution and the singular multi-pass undershooting 
solution. 

So there is the regular Vp ~ instanton, with starting 
point (j)o , and on either side of it there are singular solu- 
tions with starting points 0o i Scj). The regular instanton 
arrives at the false vacuum at t = oo. The overshooting 
solution arrives at the false vacuum too soon. However, 
the time when it does goes to infinity as 6(1) — > 0. Simi- 
larly, the undershooting solution reaches its turning point 
{(j) = 0) at a time which goes to infinity as Scf) 0. 

Now let us ask what happens to the singular solutions 



if we increase Vp while leaving the starting point fixed. 
This is a perturbation of the potential rather than of the 
initial value. If Vp is very small, then the undershoot- 
ing solution will still undershoot and the overshooting 
solution will still overshoot. By the results of Sec. IIII Dl 
between these two singular solutions must lie at least one 
regular compact instanton. 

Having demonstrated the existence of an instanton, we 
will now argue that the associated rate of decay changes 
continuously as Vp is increased away from zero. We have 
noted that the time at which the singular solutions at 
Vp = hit the singularity can be made arbitrarily large 
by choosing the perturbation dcj) to be small. By continu- 
ity of the singular solutions and by interpolation, this im- 
plies that a regular instanton exists at positive Vp whose 
size increases without bound as l^p — > 0. Continuity and 
interpolation also imply that the size of the true vacuum 
bubble is roughly the same both for the singular and 
regular solutions, and both at Vp = and infinitesimal 
positive Vp. Hence the Vp > singular solutions will 
spend all but a fixed amount of their time near the false 
vacuum^ before over- or undershooting. Therefore, the 
regular Vp > instanton will be virtually identical to the 
background de Sitter space over a volume that diverges 
as Vp 0, differing only in a "bubble volume" that re- 
mains finite as Vp — > 0. Hence, we expect the Vp > 
decay rate to be comparable to the Vp = decay rate. 

In the case where flat space is stable, continuity simply 
means that the decay rate should vanish as — > 0. 
The assumption of stability requires that any compact 
regular solution at some infinitesimal Vp > is deformed 
into a regular solution that is still compact, with finite 
action, at Vp = 0. Meanwhile, the background instanton 
action diverges as Vp — > 0. By Eq. ©, the decay rate 
approaches zero. 



IV. THE CDL SOLUTION SPACE 
A. Classification and diagrams 

In this section, we use the results derived above to 
develop a full picture of the Coleman-De Luccia solution 
space. We will be interested in how regular solutions 
can appear or disappear as the false vacuum energy Vp 
is varied beyond the infinitesimal neighborhood of Vp = 
0. We continue to distinguish the two important cases 
discussed above: whether tunneling is allowed at Vp = 
or not. 

For each of these cases, we have picked a simple poten- 
tial Vo . Figures|niand[21show solutions as a function of the 
starting point (/)o and the parameter Vp.^ More compli- 



One might imagine a situation where perturbing the starting 
point in either direction has the same effect, say undershoot. 
This is conceivable, but it is requires the first-order change to 
accidentally be zero at the instanton, which is nongeneric. 



* In fact, exponentially close. 

® This diagrammatic technique was explained to us by Matthew 
Kleban; see Batra and Kleban, to appear. 
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cated diagrams arc possible for different potentials. Our 
goal here is to analyze what appear to be the two simplest 
situations. Their structure is surprisingly rich. Here we 
will analytically explain the features they show; in the 
following section, we will confirm them numerically. 

Recall that the top of the barrier (the minimum of the 
Euclidean potential U) has been chosen to reside at </) = 
0. As long as U{0) < 0, a generic starting point {Vp, (j)o) 
will yield a singular solution. We classify each starting 
point by how many times the field passes through the 
local minimum of the Euclidean potential before escaping 
to infinity, the number of passes. These two-dimensional 
sets of points will be separated by one-dimensional lines 
corresponding to the regular solutions. 




(a) 








FIG. 2: A diagram showing solutions with starting point cl>{0) 
at the regular pole, in a potential shifted vertically by Vf. The 
numbers indicate the number of passes for each solution before 
the field escapes to infinity. The lines divide regions with 
different numbers of passes and represent regular solutions. In 
this example, tunneling is allowed for Vf = 0. The instanton 
of interest is the thick line. For positive Vf, the instanton is 
compact, and the part of the line at negative <^(0) represents 
the same solution, thinking of the opposite pole as the starting 
point. For Vf < 0, the instanton becomes noncompact so 
there is only one origin of polar coordinates. 




(b) 

FIG. 1: The Euclidean potentials we are interested range from 
that shown in (a) to that in (b). 

Note that in these diagrams a regular compact solution 
with an odd number of passes is represented by two points 
in the diagram, because either pole could be considered 
the starting point. All other solutions, in particular all 

singular solutions, are represented by one point, the value 
of the field at the regular pole. 

The fact that we can define the number of passes on a 
space consisting of the starting point and Vp , as sketched 
in the diagrams, immediately tells us that regular solu- 








FIG. 3: Here, tunneling is forbidden for Vf < 0. The single- 
pass instanton is the thick line; for a given Vf, there are two 
values of (f){0) which represent the field values at the two poles. 
For negative Vp the instanton remains compact, so it no longer 
mediates decay. 
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tions cannot simply disappear as we dial Vp; they must 
annihilate with other regular solutions in a consistent 
way. For example, a fairly general phenomenon is the 
following. As Vp is increased, raising the potential up 
to higher cosmological constant, some instantons disap- 
pear, annihilating at = 0. What is happening is that 
near the bottom of the Euclidean potential well there is 
a characteristic frequency of oscillation which is indepen- 
dent of Vp. There is also a characteristic de Sitter time 
given by the size of the instanton with the field sitting 
at the bottom- the Hawking-Moss instanton. Increasing 
Vp decreases the de Sitter time, so there is no longer 
enough time to have as many oscillations. In the figure, 
one can see the regions with multiple oscillations disap- 
pearing one by one as Vf grows until eventually only 
the Hawking-Moss instanton is left (the vertical line at 
= 0). 

In our figures, the numbers increase monotonically to- 
wards = for fixed Vp , meaning that starting points 
closer to the minimum allow more oscillations. This is a 
common situation, but as described by Hackworth and 
Weinberg 53], more interesting possibilities are not dif- 
ficult to arrange. 

Our method is extremely useful in constraining the 
presence and location of regular solutions, and it gener- 
alizes to a variety of situations. For example, one could 
consider a bigger family of deformations rather than sim- 
ply shifts of the potential. In other words, one could 
include any number of parameters on the same foot- 
ing as Vf in the analysis. In the region of parameter 
space where regular solutions are nongeneric, our method 
should apply. We consider a space in which one axis sets 
the initial condition 0(0) and the other axes represent 
the values of the parameters. The number of passes is 
a well-defined function on this space. As above, regular 
solutions generically appear at the boundaries where the 
number of passes changes. 



B. Differences between the diagrams 

Now let us discuss the differences between figures 12 
and 121 When tunneling is forbidden at Vp = (figure 
13), nothing dramatic happens to the regular single-pass 
instanton as Vf 0. That it remains compact may be 
seen by noting that at Vp = there are two regular single 
pass solution points, one at positive 4> and one at negative 
(j). These two points are really the same solution; either 
pole of the sphere can be considered the starting point. 
For Vp < our regular solutions persist, but since they 
are compact they do not describe tunneling. 

In contrast, figure |21 describes a situation where tun- 
neling is allowed at ^p = 0. This figure may look more 
fine-tuned than the one which does not allow tunneling, 
because various lines meet right at = 0. However, this 
behavior is required, based on the results of the previous 
section showing that the regular instanton is getting very 
large as Vp ^ 0. The false vacuum pole of the one-pass 



instanton (on the left of the diagram) gets closer and 
closer to the extremum. This allows the field to remain 
near the false vacuum for a very long time, resulting in a 
big instanton that looks mostly like the background false 
vacuum de Sitter solution. (This behavior results in the 
singular one-pass region on the left side of the diagram 
getting squeezed away and disappearing in the limit.) 

Also, on the right side of figure[21 the regular one- and 
two-pass solutions must approach each other as Vp — )■ 
0, squeezing the singular two-pass region. The two-pass 
instanton jSjf is a Euclidean solution which is symmetric 
about the equator, so = at the equator. For very 
large instantons, two-pass solutions require (j) to be zero 
at the equator, while one-pass solutions require to be 
miniscule at the equator so that it will be zero at the 
opposite pole. As Vp ^ and the instanton gets very 
large, the distinction between these two conditions goes 
away and the instantons merge, becoming the single-pass 
noncompact instanton for Vp < 0- 

Our diagrammatic method demonstrates other surpris- 
ing properties of the Euclidean solutions. For example, 
when tunneling is allowed, then a solution which starts 
extremely close to the false vacuum will pass the origin 
twice before escaping to infinity, while if tunneling is not 
allowed it will only pass the origin once, and escape to 
the opposite side. There is no obvious reason this prop- 
erty is related to tunneling, but our diagram suggests it 
is the case and we have verified this numerically for some 
potentials. 



C. Starting close to <^ = 

Let us explain another notable feature of our figures: 
The number of oscillations at small amplitude decreases 
as the cosmological constant increases. 

For any value of the parameters, there is a trivial solu- 
tion where the field just sits at the local minimum of the 
Euclidean potential, the Hawking-Moss instanton. Since 
we assume that this point has positive vacuum energy 
(C/(0) < 0), the geometry is a four-sphere. For suffi- 
ciently small perturbations about this instanton, we can 
compute analytically the number of passes. Because the 
field probes only the immediate neighborhood of its min- 
imum, the Euclidean potential can be approximated as a 
simple harmonic oscillator in this limit. Also, the geome- 
try can be taken to be just the Hawking-Moss instanton, 
since corrections to the geometry due to oscillations of 
the field are quadratic in the amplitude of the oscilla- 
tions. 

The Hawking-Moss instanton is 

ds^ = dt'^ + p{tfdnl (25) 

with 

p{t) = sin(i/oO . (26) 

-no 
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Here Hq is the Hubble constant determined by the value 
of the potential at the Hawking-Moss instanton, 



D. Behavior as the maximum of the potential 
approaches V — 



^(0) 



(27) 



Small oscillations are then governed by the equation 

+ 3i?o cot{Hot)^ = ~U"{Q)(I) . (28) 

Defining the dimcnsionlcss variable t — H^t, the equa- 
tion takes the form of the eigenvalue equation for the 
Laplacian on a 4-sphere,^° 



+ 3COt(f)0 : 



U"(Q) 



tt2 
^0 



(29) 



where (f> is now the derivative with respect to t. The 
quantity U"{0)/Hq plays the role of the eigenvalue. 

Regular solutions exist only if the eigenvalue is of the 
form n(n + 3) for a positive integer n. These solutions 
are spherical harmonics which are homogeneous on the 
three-sphere. They have n zeroes, so for our purposes 
they are n-pass solutions. Of course the approximation 
we are using is only valid for infinitesimal oscillations, 
so the physical statement is that as we shift the poten- 
tial, infinitesimal-amplitude regular solutions only exist 
at special values. This is consistent with our general rule 
that regular solutions are nongeneric. 

The number of passes will change as we shift the po- 
tential; Hq changes if we add a constant Vp to the po- 
tential, while C/"(0) is unaffected. Recall that we denote 
the vertical shift of the potential by Vp, so even though 
these solutions do not care about the value of the en- 
ergy in the false vacuum, in our notation we say that 
infinitesimal-amplitude solutions exist at special values 
of Vp. Between values of Vp allowing a regular solution 
with n passes and a regular solution with n + 1 passes, 
the singular solutions have n -I- 1 passes. So for singular 
solutions the number of passes is given by the smallest 
integer n such that 



n{n + 3) > 



U"{0) 

^0 



(30) 



For sufficiently large positive Vp, it is clear from the 
formula that the number of passes will be n = 1. This 
happens because the Hubble time for the Hawking-Moss 
instanton becomes short compared to the period of the 
harmonic oscillator. 

Once again, we see that regular solutions exist at the 
boundary between singular solutions with different num- 
bers of passes. As explained by Gratton and Turok [H^j ]. 
the number of passes is equal to the number of negative 
modes for perturbations around the Hawking-Moss solu- 
tion. (The number of negative modes is important for the 
interpretation of the instanton, but in this paper we are 
mainly interested in the structure of the solution space.) 



So far, we have assumed that the maximum of the 
potential (the Hawking-Moss point) has positive energy: 
V{0) > 0. What happens if we allow ¥(0) = 0? This 
corresponds to the local minimum of the Euclidean po- 
tential moving up to ?7 = 0. The results are dramatic. 
We can no longer argue that solutions are generically 
compact, or generically singular, because there is a new 
possible asymptotic behavior: the field can asymptoti- 
cally approach = while the geometry continues to 
grow (p oo). This is possible because the potential 
is now zero at the minimum, so asymptotically flat solu- 
tions are possible. 

Furthermore, this asymptotic behavior is stable under 
small perturbations, since the field approaches a local 
minimum of the Euclidean potential. So generic initial 
conditions can result in noncompact, regular solution. 
Specifically, all starting points with Euclidean potential 
energy less than Up , the Euclidean potential of the false 
vacuum, lead to this asymptotically flat behavior; de- 
pending on parameters, a somewhat bigger set of start- 
ing points leads to asymptotically flat solutions. In fact, 
for the potentials we will consider numerically any start- 
ing point between (f>p and (pT will lead to a noncompact 
solution with (jj —^ 0. 

To see that this is reasonable, recall the FRW equation 



Stt 
3A/2j 



(31) 



Once again, we thank Matt Kleban for pointing this out. 



As long as the field does not escape the region of inter- 
est, H cannot go to zero except at p = oo because the 
Euclidean potential U is nonnegative in this region. If 
H cannot go to zero then it cannot become negative, so 
no anti-friction is available. So solutions starting with 
energy less than Up have no hope of escaping because no 
anti-friction is available until after they escape, and they 
would need to gain energy in order to escape. Solutions 
beginning near the true vacuum, with more potential en- 
ergy than [/p, can hope to overshoot on the first pass, 
but whether this is possible depends on parameters. 

There must be a dramatic signal of this lurking non- 
compactness as the maximum of the potential approaches 
zero from above. One place we can see this is in our 
formula for small oscillations around the Hawking-Moss 
solution. Note that in this limit the Hubble constant Hq 
of the Hawking-Moss instanton approaches zero because 
the cosmological constant of the Hawking-Moss solution 
is approaching zero. According to our formula Eq. H30|) 
the number of passes approaches infinity as Hq — s- 0. The 
reason is that as Hq the Hubble time goes to infinity 
while the period of the harmonic oscillator stays fixed, so 
there is time for an infinite number of oscillations. 

We have thus discovered an accumulation point, with 
an infinite number of regular solutions appearing at 
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Hq 0. This is shown schematically in figures |21 and El 
and will also appear in our numerical data. 



E. Starting close to the vacua 

Here we explain the behavior seen in both Fig. [21 and 
Fig. [31 for starting points very close to the false vacuum. 
We will also be able to explain some aspects of the be- 
havior just to the left of the regular single-pass instanton 
and just to the right of the true vacuum. Note in the 
diagrams that the number of passes is one for solutions 
whose starting point is just to the right of the false vac- 
uum as long as Vp > 0. However, as Vp becomes more 
negative the number of passes increases. 

It is clear that the number of passes must be one for Vp 
positive. The reason is that there is a regular compact 
solution with the field sitting on the false vacuum- just 
the usual false vacuum de Sitter solution. Also, starting 
points to the left of the false vacuum must lead to zero 
passes by our assumptions about the potential. Since the 
number of passes must change by one across a compact 
regular solution, starting points just to the right of the 
false vacuum must lead to one-pass solutions. 

This logic no longer holds once Vp < 0. There is still 
a regular solution where the field sits on top of the false 
vacuum forever, represented by the vertical line in the 
diagram, but this solution is now noncompact. We have 
no general rules about how the number of passes can 
change across a noncompact regular solution. 

To see that starting points just to the right of the false 
vacuum might result in a large number of passes for Vp 
negative, imagine a situation where Vp is very negative 
so that both vacua are at very negative cosmological con- 
stant. (Note that we are no longer interested in the limit 
of the previous section where the Hawking-Moss solution 
becomes noncompact; here we assume that it is com- 
pact.) If the field starts very close to the false vacuum, it 
will sit there a long time before beginning to oscillate, so 
that it essentially begins to oscillate at large p compared 
to all other scales. 

Recalling equation H14|) . the change in the Hubble pa- 
rameter is given by 



(32) 



We ignore the second term because we are at large p. 
The change in H over one oscillation, given by the above 
formula, is determined by the potential. It is basically 
constant as we shift the potential vertically. Now before 
the field starts oscillating the Hubble parameter is posi- 
tive, and its magnitude is determined by the cosmological 
constant at the false vacuum, so its magnitude is large for 
Vp ^ 0. As a result, one oscillation of the field produces 
a small change in AH/H. 

We have shown in section UlI Al that any solution, sin- 
gular or nonsingular, starting near the false vacuum must 
be compact. Since the solution must be compact, H 



must become negative. But since one oscillation leads 
to a small change in AH/H, it will take many oscilla- 
tions before H can become negative. Further, the field 
experiences no anti-friction as long as > 0, so it is 
guaranteed to keep oscillating until H becomes negative. 

Once H becomes negative, more oscillations are nec- 
essary before the field can escape (recall that generically 
the field escapes to ±oo), because the magnitude of H is 
bounded below by the vacuum energy. 



H^ > 



Stt 



u{4>) 



(33) 



Another way to think about it is that the field has lost 
energy to Hubble friction during its oscillations, and it 
must undergo more oscillations during the phase of anti- 
friction to recover enough energy to approach one of the 
vacua. These two points of view are related by the FRW 
equation in the limit p ^ oo. 



H' 



Stt 



E 



(34) 



where E — ^cjp + U{(j>) is the Euclidean energy of the 
field. 

Generically, the field will not return precisely to one of 
the vacua after these oscillations; it will escape with the 
radius is still extremely large. 

To summarize, for starting points just to the right of 
the false vacuum, the generic solution consists of a very 
large region where the field is essentially in the false vac- 
uum. When the field begins to oscillate there may be a 
large number of oscillations at large p, after which the 
field escapes and the radius returns to zero. The field 
generically escapes at large radius. 

Since there is a separation in scales between the size of 
the false vacuum region and the characteristic period of 
oscillation, these solutions should be well captured by the 
thin wall approximation. The thin wall approximation 
will allow us to make more concrete statements about 
the number of passes. In order to capture the poten- 
tially large number of oscillations, we slightly generalize 
the thin wall approximation: we allow ourselves to stack 
multiple domain walls on top of each other, each domain 
wall representing one pass. 

We want to construct, within the thin wall approxi- 
mation, a solution (singular or nonsingular) which has 
the features described above. It should have a regular 
pole surrounded by an enormous false vacuum region sur- 
rounded by a stack of domain walls at infinity. On the 
other side of the stack of domain walls, the radius should 
decrease back to zero. We want to allow solutions with 
a singularity at one pole as usual. Within the thin wall 



We emphasize that nowhere else in this paper do we appeal to a 
thin wall limit. 



12 



approximation, the only possible singularity is a conical 
deficit. With a conical deficit, the metric becomes 

ds^ = dt^ + fp{tfdnl , (35) 

where p{t) takes the usual Euclidean AdS form, p{t) = 
Rsmh{t/ R), and / > 0. (For / = 1 there is no conical 
deficit or excess.) 

To summarize, we are seeking a thin wall solution 
which has a regular pole surrounded by an enormous false 
vacuum region. The false vacuum region is surrounded 
by a stack of domain walls. On the other side of the do- 
main walls, we can have a region of true or false vacuum 
which generically has a conical deficit. We expect that 
the thin-wall solution is a good representation of the full 
solution in the false vacuum; we also expect it to capture 
correctly the number of field oscillations required before 
the field has enough energy to escape. However, as men- 
tioned above the field generically escapes at large radius, 
and after this point it is not near either of the vacua, so 
the thin wall approximation must be wrong. Neverthe- 
less, the thin wall approximation should be a good one 
for predicting the number of oscillations. 

The Israel junction condition relates the jump in ex- 
trinsic curvature across the stack of domain walls to the 
tension. For the situation we are considering, in which 
the radius decreases as we move away from the domain 
wall in either direction, it is 

where i?i and i?2 are the AdS radii on either side of the 
wall, /i and /2 are related to the conical deficits, a is 
proportional to the tension of one domain wall, and n is 
the number of domain walls in the stack. The left side of 
the equation is bounded below by 1/ R1 + I/R2 ■ Since the 
first pole is nonsingular by construction, we set /i = 1. 

1. Vanishing false vacuum, energy 

We first consider the special case where the false vac- 
uum has zero cosmological constant, i?F = 00. Above, 
we claimed that if the false vacuum has very negative 
cosmological constant, the number of passes (the num- 
ber of domain walls in the thin wall approximation) will 
be large. On the other hand, if the false vacuum has zero 
cosmological constant we will see that one or two domain 
walls is sufficient. The simplest solution would be a large 
region of false vacuum, one domain wall at large radius, 
and a region of true vacuum on the other side. In this 
case the left side of Eq. H36|l is bounded below by l/i?T. 
If the tension is big enough, a > 1/Rt, a solution exists. 
In this case, for sufficiently large p, we can always solve 
the junction condition for the conical deficit parameter 
/• 

On the other hand, if cr < 1/Rt, no solution exists. In 
this case, we will need two domain walls to solve the junc- 
tion condition. With an even number of domain walls we 



have false vacuum on both sides. The junction condition 
becomes 

- + ^ = 2a (37) 
P IP 

This equation can be solved for any tension cr, because 
the left side can take any value. For fixed sufficiently 
large p we can again solve for /. 

So we need one domain wall if cr > 1/Rt and two 
domain walls if cr < 1/ Rt But cr < 1/ R^ is precisely the 
condition that the false vacuum can decay! Equivalently, 
it is the condition that a noncompact regular instanton 
exists. The conclusion is that if tunneling is forbidden 
(Fig. O, starting points just to the right of the false 
vacuum have one pass, while if tunneling is allowed (Fig. 
|2Jl they have two passes. 

Another way to understand that two domain walls are 
required when tunneling is allowed is to note that Eu- 
clidean AdS space has finite extrinsic curvature at infin- 
ity. Hence, joining an AdS region to a fiat region across 
a three-sphere requires a jump in the extrinsic curvature 
that remains finite even as the location of the domain wall 
p goes to infinity. If tunneling is allowed, the tension is 
too small to account for this finite change. 

Now let us move to the right part of Fig. [3 We can in- 
fer the behavior of solutions which start just to the left of 
the single-pass noncompact instanton, that is, solutions 
which barely undershoot the false vacuum. The regular 
single-pass instanton achieves the false vacuum at very 
large p; a slightly undershooting solution will spend a 
long time close to the false vacuum before falling back. 
Since we have just analyzed solutions which spend a long 
time and grow to large radius near the false vacuum, the 
subsequent evolution will be the same. 

In the thin wall approximation, the equivalent state- 
ment is that the single-pass instanton and the pure false 
vacuum solution have the same asymptotics, so we can 
replace the false vacuum region in our previous solution 
by the single-pass regular instanton. The conclusion is 
that a solution starting just to the left of the regular 
instanton will have one more pass than a solution which 
starts just to the right of the false vacuum. This explains 
why the number of passes jumps from one to three across 
the regular instanton. 



2. Negative false vacuum energy 

Next, let us consider the lower part of Fig.[21and Fig.O 
the region with Vp < 0. Again we folcus on starting 
points just to the right of the false vacuum. Now it is no 
longer clear that we can have a solution with two domain 
walls, because even the false vacuum has finite extrinsic 
curvature as p ^ 00; as we argued below Eq. 1)321) a 
large number of domain walls may be necessary. Two 
types of solutions are possible: with an even number of 
domain walls, we have false vacuum on both sides, and 
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for a solution to exist we need 



l/i?F + l/i?T < na, n even. 



(38) 



For an odd number of domain walls, we have true vacuum 
on one side, and for a solution to exist we need 



l/i?F + l/i?T < ncr, 



odd. 



(39) 



The minimum number of passes is determined by the 
smallest number n of domain walls such that one of the 
inequalities is satisfied. It is clear that as Vpbecomes 
more negative the required number of passes n becomes 
larger. This is reflected in both diagrams. 

Additionally, if tunneling is allowed for a given < 0, 
then the minimum number of passes always occurs for n 
even. This is because the condition that tunneling is 
allowed is I/Rt — i/Rp < c. When this is satisfied, 
equation H38|) is always satisfied for smaller n than equa- 
tion l|39|l . The result is that along the left side of the 
diagram, all numbers should be even for Vf < as long 
as tunneling is allowed. 

Let us move to the right of the figures. As in the Vp = 
case, we parlay our results about solutions starting near 
the false vacuum into information about starting points 
just to the left of the single-pass regular instanton. Once 
again the number of passes just to the left of the noncom- 
pact instanton should be one greater than the number of 
passes for a starting point just to the right of the false 
vacuum for fixed Vp. 

Hence, if tunneling is allowed, only odd numbers 
should appear immediately to the left of the regular in- 
stanton. On the other hand, if tunneling is forbidden 
then as we dial Vp the number of passes n near the false 
vacuum will change by one at a time, though there will be 
much more parameter space where the number is even. 



V. NUMERICAL EVIDENCE 
A. Method 

In order to exhibit the features discussed in the previ- 
ous sections, we turn to numerical solutions of the CDL 
instanton equations l(2Jl. 

Following Ref. J£J, we rescale the physical quantities 
to yield dimensionless variables. The potential for (jj can 
be written 



where 



vi4>) = ii\{x) 



X = (f)/M 



(40) 



(41) 



Here, M is the scale over which the potential has non- 
trivial features, and fx^ is a characteristic energy density. 
We assume that <C M < Mpi where Mpi is the Planck 
scale. The two mass scales have been extracted so that v 




FIG. 4: A graph of the potential v{x) for & = 1, and z 
{1,0,-1}. 
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FIG. 5: A graph of the potential v{x) for b = 0.3, and z = 
{1,0,-1}. 



is a function that can be approximated as a polynomial 
with coefficients of order one. 

We employ the same quartic potential used in Ref. , 
which takes the form 



v{x)^f{x)-{l + z)f{xp) 



where 



X 

T 



bx^ 



X 

T 



(42) 



(43) 



The function f{x) has two negative local minima at 
and XT, with f{x-p) > /(xt), and a local maximum at 
X = 0. The potential v(x) share these properties. 

The adjustable parameter z allows us to tune the 
false vacuum energy by shifting the entire potential with 
v{xf) oc z. Thus, z plays the role of the shift Vp in 
the previous sections. In particular, the interesting limit 
Vp corresponds to z ^ 0. (Note that /{xp) < 0, so 
z > corresponds to Vp > 0.) 

The parameter b, which we take to be strictly positive, 
controls both the width xt — xp and relative heights of 
the vacua v{xf) — v{xt)- Figures ^ and [S] show plots of 
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v{x) of b = 1 and 0.3, respectively, for various values of 
z. 

The radius and time can also be made dimensionless 
by rescaling by appropriate powers of the mass scales fi 
and M: 



M ' 

Id ' 
M 



(44) 
(45) 
(46) 



We have also defined a dimensionless quantity e control- 
ling the strength of gravity. For e <g; 1, gravity has a 
negligible effect on the decay rate, but when e is of or- 
der one, gravity can be important. For example, it can 
completely suppress the decay of flat space. 

Let us rewrite the Euclidean CDL equations Q and 
© in terms of these dimensionless variables: 



X + 3—x = v' (x) , 
r 



1 2 2 

1 = e r 



X 

Y 



v{x) 



(47) 
(48) 

(49) 



where an overdot (prime) denotes differentiation with re- 
spect to s {x). 

Using Mathematica, we numerically integrate these 
equations with initial boundary conditions 



r(0) ^ 1 
i(0) ^ 



r(0) = , 
a;(0) xq 



(50) 
(51) 



to yield functions r(s) and x{s). The initial position xq 
along with z, b, and e form a set of four adjustable pa- 
rameters on which the solutions depend. 

By numerically constructing many solutions, both sin- 
gular and regular, for a variety of parameters, we will 
be able to verify the assertions made in sections IlIII and 
IIVI For generic choices of parameters, the numerical 
integration produces compact singular solutions where 
the field x{s) inevitably escapes to ±00, as explained in 
Sec. nil CI For regular compact solutions, the unbounded 
anti-friction as r — s- makes the numerics difficult to con- 
trol near the far pole. By tuning xq, however, regular 
solutions can be very well approximated. The range of 
r over which x{s) can be made to remain near xt or x-p 
before rolling off to ±00 is limited only by the tuning of 
Xq and calculational precision. 

Noncompact solutions can likewise be approximated by 
fine-tuning xq- In accordance with the arguments of sec- 
tion [ffll these are necessarily regular one-pass solutions 
starting at the true vacuum and asymptoting to the false 
vacuum. As with the approximately regular compact so- 
lutions, the inevitable result of imperfect tuning of xq 
is that x{s) eventually either over- or under-shoots x-p, 




7.5 10 12.5 15 



FIG. 6: This graph shows a family of one-pass solutions with 
b = 0.3, e = 0.6, and z = {1,0.1,0.01,0}. As z 0, these 
compact solutions smoothly approach the z — solution, 
which for this value of (6, e) is also compact. The decay of 
flat space is forbidden. This reproduces the results of 



leading to a singularity. However, with sufficient tuning, 
we can obtain a large (in s) region where r = 1 and reach 
a large maximum r before the singularity. 



B. The decay of nearly flat space 

Our first goal is to confirm the result of Sec. IIIII that 
the decay rate is continuous as 2: ^ 0. 

As was pointed out in Ref. 48], for certain potentials 
a seemingly metastable Minkowski space (z = 0) can be 
stabilized by gravitational effects. To exhibit this effect, 
we tune e which controls the strength of gravitational 
corrections. Keeping 6 = 1 fixed, we find a critical value 
of e « .35. Below this value, the z — one-pass solu- 
tion is noncompact, implying the instanton corresponds 
to an actual decay of Minkowski space. Above the crit- 
ical value, the regular one-pass solution is compact at 
z = and therefore does not represent decay. A simi- 
lar transition occurs when varying b for fixed e. For fixed 
e = 0.6, for example, the z — regular one-pass solutions 
are compact below b « 0.74 and noncompact above this 
value. 

We now turn to the question discussed in Sec. IIIII the 
behavior of compact, regular one-pass z > solutions in 
the limit where z — s- 0. The asserted continuity of the 
decay rate translates into different behavior, depending 
on whether the false vacuuum at z = is stable or not. 
We chose e = 0.6, & = .3 to study the stable case and 
e = 0.6, 6 = 1 for the unstable case. 

For each case, we probed the z —^ limit by numeri- 
cally computing a family of regular single-pass solutions 
with z — {1,0.1,0.01,0}. The stable case is shown in 
figure El As expected, the family of compact solutions 
smoothly approaches the compact z = solution. This 
is the behavior anticipated in Sec. IIIII 
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In the unstable case, the compact z > solutions dra- 
matically grow in size as z — > 0, as shown in figure \7\ 
For z < 0.1 we lack the numerical precision to follow the 
solution all the way to the far pole at s > 30 because 
the corresponding dS spheres are growing so large. For 
z = 0.1 at least the equator where r = 0ats = 22is 
still within our computational range. We can clearly in- 
fer, though, that as z ^ the compact instantons are 
growing so as to reach infinite size in the limit. The 
z — 0.01 and z = solutions are virtually indistinguish- 
able. Again, this supports the behavior anticipated in 
Sec.lml 

In summary, our analytic argument that the decay rate 
is continuous is borne out by the numerical evidence. 



C. Plots of the solution space 

Our next numerical goal is to verify the broader anal- 
ysis of the CDL solution space provided in section IIVI 
For fixed e — 0.6, we considered the two important cases 
b — 1 (flat space decays) and 0.3 (flat space is stable). 
For each case, we computed the number of passes for 
10^ solutions with z G [—1,9] and xq e [xp,xt]- The 
resolution is approximately Sxq ~ 10~^ and Sz ^ 10^^. 

The results are shown in flgures|Sland|51 The number 
of passes for each region is labeled. The curves demarcat- 
ing the boundaries between regions with different num- 
bers of passes are the locations of the regular solutions. 

To better access features near xp and xt, we also ran 
10^ points close to each edge for each value of b, with 
exponentially decreasing step size in xq as the edge was 
approached. Our resolution was then 5\og{xo) ~ 10~^ 
and Sz ^ 10^^. The magnified left edges of figures |H1 and 
Oare shown in figures ITUlandlT^ Likewise, figures [TTland 
1131 zoom in on the right edge. 



FIG. 8: The {xo,z) solution space for 6=1 and e = 0.6. 
This numerical plot reproduces the general features of figure 
121 Here, tunneling is allowed at z = 0. The numbers to 7 
denote the number of passes for solutions in a given region. 
The curves at the boundaries between regions are regular so- 
lutions. The vertical lines are the trivial solutions at xt and 
xp, and the thick blue line is the single-pass instanton. The 
Hawking-Moss solution (a vertical line at a; = 0) is hidden by 
the axis. Features near xt and xf are not well-resolved; see 
figures El and fTTI 



D. Discussion 

A number of important features predicted in figures 
|51 and 13 are reproduced numerically in figures |H1 and |51 
As we argued in Sec. IIII Dl when crossing any regular 
compact solution, the number of passes should jump by 
one. This can be seen in figures |H1 and [HI It can also 
be seen that curves representing noncompact solutions 
separate singular regions whose number of passes differ 
by as much as seven. 

Another notable feature of the z ^ limit in the un- 
stable case (figure |2) was the merging of solutions as 
z —> 0. As explained in section IIV Bl the regular one- 
and two-pass solutions on the right merge into each other 
at z = 0, as do the zero- and one-pass solutions on the 
left. This was the most obvious feature distinguishing 
the unstable from the stable diagram (figure |21). 

This merger can be roughly made out in our numerical 
plot for the unstable case, figure |S1 and it is notably 
absent in the stable plot, figure El This is best seen in 
the refined plots for the left and right edge; note how the 
thick blue line passes through z = in figures 1121 and 
[T^ . This supports our arguments in Sec. IIV Bl 

However, neither figure O nor its refinements, figures 
[TUl or ^2 are able to show that the merger takes place 
exactly at z = 0. The curves become so close together 
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FIG. 9: The (1:0,2) solution space for b — 0.3 and e — 0.6. 
This numerical plot reproduces the general features of figure|3 
Tunneling is not allowed at z = 0. The numbers to 7 denote 
the number of passes for solutions in a given region, and the 
curves represent both regular solutions and the boundaries 
between regions. The vertical lines are the trivial solutions at 
XT and xf, and the thick blue line is the single-pass instanton. 
Features near xt and xp are not well-resolved; see figures fT^ 
and [131 



FIG. 11: A close-up of the lower right corner of fig. m (t un- 
neling allowed ai z = 0) showing the numerically computed 
{xo,z) solution space with 6=1 and e = 0.6. The regular 
one- (thick, blue line) and two-pass solutions merge as z de- 
creases. Due to our limited precision, they appear to meet at 
z = 0.3 rather than at z = as expected. 
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FIG. 10: A close-up of the lower-left corner fig. H (tunnel- 
ing allowed at 2 = 0), showing the numerically computed 
(xq, z) solution space with b = 1 and e = 0.6. The thick, blue 
one-pass curve is nearly vertical and is merging with the xy 
solution. Although we expect that they meet at 2 = 0, our 
limited precision makes them appear to merge at higher z. 



FIG. 12: A close-up of the lower left corner fig. El (tunneling 
forbidden at 2 = 0) showing the numerically computed [xq, z) 
solution space with 6 = 0.3 and e = 0.6. One can see that 
the one-pass regular solution passes through 2 = without 
merging. 
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FIG. 13: A close-up of the lower right corner fig. |2l (tunneling 
forbidden at z = 0) showing the numerically computed (so, z) 
solution space with h = 0.3 and e = 0.6. Again, one can 
see that the one-pass regular solution passes through z = 
without merging. 



For example, in the 6=1 case, at z = 1 the right-hand 
side of Eq. H52|l is 18.3 which yields n = 3. From figure 
IHl we can see the curve of regular three-pass solutions 
crosses the origin just above z — \. Below it, in particular 
at z = 1 near x — Q, are singular three-pass solutions. 

Confirming a prediction from Sec. IIVDI the regular 
solution curves begin to pile up when z — > —1, in both 
figure IHl and figure El Ultimately, the resolution is insuf- 
ficient to detect every transition between constant-pass 
regions. As z — > — 1 the number of oscillations increases 
rapidly, the amplitude decreases, and the size of the so- 
lutions grow very large, s ^ 1. For example, when 6=1 
, z — 0.99, and \xq — xtI = 10~^, x{s) oscillates for so 
long around the Hawking-Moss solution, that the numer- 
ical precision is insufficient to follow the solution back to 
r = 0, and instead errors build up and lead to a singu- 
larity at s ~ 130. 

The plots also connect with the discussion in Sec. IIVI 
For z < 0, the regular solutions appear to merge in pairs 
before reaching x-p or xt. However, we can resolve the 
two curves, and, when we zoom in on the edges by plot- 
ting figures IIUI and 1121 logarithmically in xq — in or 
figures 111! and 1131 logarithmically in xt — a;o , we can see 
the curves reach xp or xt at separate points. 



that even at z = 0.1, a resolution of 6xq ^ 10^^° is 
insufficient to detect the region between them. 

The plots also confirm the estimate in Sec. IIV CI for 
the number of oscillations for solutions near xq = 0: The 
regular solution with the largest number of passes at a 
given z > — 1 is given by the largest integer n satisfying 
Eq. H3U|) . which in dimensionless variables is 
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